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1.  INTRODUCTION 


The  gases  discharged  from  a  gun  muzzle  after  firing  determine  the 
strength  of  the  blast  which  affects  firing  crew  safety  and  shell  accuracy. 

In  order  to  understand  this  phenomenon,  a  number  of  theoretical  and  experi¬ 
mental  studies  have  been  conducted.  Extensive  investigations  of  muzzle 

1  2 

blast  were  carried  out  by  Schmidt  and  his  associates  at  BRL.  A  highly 
complex  time -dependent  flow  field  is  noted  in  their  flow  visualizations  results. 
For  example,  a  schematic  of  the  flow  pattern  is  depicted  in  Fig.  1.  It 
illustrates  a  large  overexpansion  jet,  recompression  shocks,  contact  surface, 
Mach  disc,  and  free  air  blast.  In  early  analytic  work,  Oswatitsch  modeled 
the  initial  unsteady  phase  of  this  type  of  flow  as  a  spherical  blast.  He  applied 
the  method  of  characteristics  to  calculate  the  blast  field  about  the  gun  muzzle. 
Subsequently,  he  noted  that  in  the  region  between  the  muzzle  and  the  inward 
facing  shock  one  can  approximately  model  the  flow  by  a  steady  jet  theory. 
Erdos  and  Del  Guidice^  have  used  this  suggestion  to  evaluate  the  muzzle  flow 
properties  along  the  symmetry  line. 


^Schmidt,  E.M.  and  Shear,  D.  D.,  "The  Flow  Field  About  the  Muzzle  of  an 
M-16  Rifle,  "  BRL  Report  No.  1692,  U.  S.  Army  Ballistic  Research  Labora¬ 
tory,  Aberdeen  Proving  Ground,  Md.,  Jan.  1974. 

^Gion,  E.J.  and  Schmidt,  E.M.,  "Pressure  Measurements  on  a  Muzzle 
Brake  Simulator,  "  J.  of  Ballistics,  Jan.  1978. 

^Oswatitsch,  K.,  "Flow  Research  to  Improve  the  Efficiency  of  Muzzle  Brakes, 
Part  I  through  Part  III.  "  Muzzle  Brakes,  E.  W.  Hammer  (ed.),  Franklin 
Institute,  1949. 

"^Erdos,  J.  and  Del  Guidice,  P.,  "Calculation  of  Muzzle  Blast  Flow-Fields,  " 
AIAA  J.,  J3,  1975,  p.  1048. 

(This  study  was  supported  by  the  U.  S.  Army  ARADCOM  under  Space  Division, 

Air  Force  Systems  Command  Contract  No.  F04701-79-C-0080. ) 
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Muzzle  Flow  Schematic 


Recently,  several  attempts^ have  been  made  to  apply  numerical 
methods  to  evaluate  the  fully  two-dimensional  or  axisymmetric  time- 
dependent  muzzle  flow  field.  Numerical  calculations  can  predict  the  muzzle 
flow  development  in  greater  detail  than  in  earlier  analytic  models.  How¬ 
ever,  none  of  the  calculations  mentioned  appear  to  have  been  compared  to  the 
existing  experimental  data  in  order  to  demonstrate  their  accuracy  and 
validity. 

As  mentioned  before,  the  muzzle  blast  flow  field  is  complex  and 
inherently  unsteady.  It  provides  a  severe  test  of  any  numerical  simulation 
technique  due  to  the  wide  spectrum  of  flow  conditions  which  exists  as  the 
transient  flow  field  develops.  However,  the  detailed  simulation  of  the  prob¬ 
lem  would  considerably  aid  attempts  to  alleviate  undesirable  weapon  charac¬ 
teristics  such  as  recoil,  noise,  flash,  and  projectile  dispersion.  Until  now, 
devices  to  suppress  these  undesirable  characteristics  are  essentially 
designed  by  trial-and-error  experimental  programs. 


^Taylor,  T.  D.,  "Calculation  of  Muzzle  Blast  Flow  Fields,"  PA-R-4155, 
Picatinny  Arsenal,  Dover,  N.J.,  Dec.  1970. 

^Moore,  G.  R.,  "Finite  Difference  Calculations  of  the  Free-Air  Gun  Blast 
About  the  Muzzle  of  a  5"/54  Naval  Gun,  •'  Naval  Weapons  Laboratory, 

Report  TR-2794,  Sept.  1972. 

^Ishiguru,  T.,  "Finite  Difference  Calculations  for  Two-Dimensional  Unsteady 
Expanding  Flows.  "  AIAA  J .,  _1^(2),  Feb.  1972. 

^Maille,  F.  H.,  "Finite  Difference  Calculations  of  the  Free-Air  Blast  Field 
About  the  Muzzle  and  a  Simple  Muzzle  Brake  of  a  105  mm  Howitzer,  " 

Naval  Weapons  Laboratory,  Report  TR-2938,  May  1973. 

^Traci,  R.M.,  Farr,  J.  L.,  and  Liu,  C.Y.,  "A  Numerical  Method  for  the 
Simulation  of  Muzzle  Gas  Flows  with  Fixed  and  Moving  Boundaries,  " 
BRL-CR-161,  U.S.  Army  Ballistic  Research  Laboratory,  Aberdeen  Proving 
Ground,  Md.,  June  1974. 

^*^Moretti,  G.,  "  Muzzle  Blast  Flow  and  Related  Problems,  "  AIAA  paper 
78- 1190',  Seattle,  Wash.,  July  1978. 
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In  this  report,  we  investigate  the  phenomenon  of  muzzle  blast  by 
extending  the  numerical  procedure  used  originally  in  Ref.  5.  In  this  pro¬ 
cedure  we  apply  Godunov's  method  along  alternating  directions  to  solve  the 
inviscid  equations.  The  procedure  works  well  except  along  the  axis  of  sym¬ 
metry  where  special  precautions  are  necessary  due  to  the  extreme  gradients 
in  the  flow.  Details  of  the  application  and  calculations  results  are  presented 
in  the  following  sections. 
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2.  FORMULATION 


In  order  to  analyze  the  flow  from,  the  muzzle  of  a  gun,  we  selected  an 
inviscid  flow  model.  Such  a  flow  is  described  by  the  axisymmetric  time- 
dependent  Euler  equations  which  can  be  written  in  the  form 


(py)t  +  (pvy)y  +  =  o 

(puy)^  +  [y(p  +P^^)]^  +  (puvy)y  =  0 
(pvy)^  +  (puvy)^  +  [y(p  +  pv^)]^  -  p  =  0 
(pEy)^  +  (puHy)^  +  ^vHy)^  =  0 


where 

E 

e 

H 

t 

X 

y 

p 

P 

u 

V 

p 


e  +  1/2  (u^  +  v^) 
specific  energy  of  the  fluid 
E  +  P/P 
time 

direction  along  the  axis  of  the  gun 

radial  direction  measured  from  the  axis  of  symmetry 

pressure 

density 

X  velocity 

y  velocity 

(y-i)pe 
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Godunovas  scheme  was  applied  to  discretize  this  gov'erning  system; 
however,  we  found  that  the  nonlinear  form  of  Godunovas  method  must  be  used, 
since  the  linear  version  of  the  scheme  was  discovered  to  be  numerically 
unstable  in  regions  of  extreme  gradients.  Also,  an  ADI  procedure  was  em¬ 
ployed.  The  split  equations  took  the  following  form: 

Set  I  (for  integration  in  the  x  direction) 


pn+l/2  _  pn  _  1 

hi  X  L  i,J  1-1, jJ 

=  (p^)l.  -  r/h^  [(P  +  Ru2)._  .  -  (P  + 

KY''"  =  Ki  ] 

=  tf>E)"  .  -  r/h^  [ruh._  .  -  J 


Set  II  (for  integration  in  the  y  direction) 


►  -  r/h  r  RV.  .  -  RV.  .  ,  1 

1,  J  i,J  yL  1,  j  1,  j-lJ 


,  .n  +  l  ,  ,n+l/2 

(nu).  .  =  tou).  .  -  r/h 

^  1,  J  ^  1,  J  > 


RUV.  .  -  RUV 
1,  J 


i,  j-1- 


RUV.  .  +  RUV. 

1,  J  1, 


T 


Godunov,  S.K.,  "Finite  Difference  Method  for  Numerical  Solutions  of  Equa¬ 
tions  of  the  Equations  of  Fluid  Dynamics,  "  Mat.  Sbornik,  47(89),  No.  3, 
1959.  ~ 
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(pE)f/  =  <^E)U  -  ’■%  ['"'"'“’i,  j  ■  j-l] 

In  these  expressions  each  grid  has  an  ordered  pair  of  subscripts  (i,  j),  with  i 
denoting  the  x  direction  and  j  the  y  direction.  The  capital  quantities  (e.g., 

R,  U,  V,  H,  and  P)  refer  to  values  of  the  flow  variables  on  the  cell  edge  which 
are  computed  by  the  nonlinear  forms  of  Godunov's  scheme  (see  Fig.  2). 

Figure  3  illustrates  the  geometry  and  boundary  conditions  used  in  our 
formulation.  In  the  present  approach  the  shell  is  included  in  the  analysis. 
Boundaries  A,  B,  and  C  are  solid  boundaries  with  zero  normal  velocity, 
while  boundary  D  is  the  shell  base  with  a  prescribed  velocity  u^.  Boundary 
E  is  an  expanding  grid  line  which  moves  outward  just  ahead  of  the  muzzle 
blast.  The  grid  is  allowed  to  expand  to  a  specified  limit  in  x  and  y;  then, 
"flow-through"  boundary  conditions  are  imposed  along  the  line.  This  amounts 
to  requiring  each  cell  outside  the  boundary  to  have  the  same  flow  conditions 
in  the  outward  normal  direction  as  the  cell  inside  the  boundary. 

Boundary  F  is  the  axisymmetric  line  across  which  there  is  no  flow, 
and  G  is  the  muzzle  exit  condition.  Initially,  the  conditions  at  G  were  speci¬ 
fied  as  the  muzzle  launch  conditions  for  all  time.  However,  this  was  found 
to  be  incompatible  with  the  real  problem,  since  extreme  pressures  were 
generated  in  the  muzzle  blast  which  eventually  caused  the  computations  to 
become  unstable.  When  the  condition  on  G  was  relaxed  and  the  actual  flow 
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Evaluation  of  Boundary  State  in  Godunov  Scheme 
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Fig.  3.  Geometry  for  Muzzle  Blast  Calculation 


within  the  gun  was  computed  or  specified  in  conjunction  with  the  muzzle  blast 
expansion,  the  results  were  very  satisfactory.  In  the  final  configuration, 
the  grid  system  considered  flow  both  inside  and  outside  the  muzzle  which  is 
designated  as  boundary  A. 

The  initial  conditions  for  the  problem  consisted  of  ambient  atmospheric 
conditions  within  the  boundary  AEC  at  launch.  Within  the  gun  ABFG,  the 
conditions  are  prescribed  as  the  launch  conditions  and  must  be  provided  by 
an  independent  interior  ballistic  calculation. 

The  velocity  of  the  boundary  D  is  prescribed  as  the  shell  launch  velocity 

u  .  This  boundary  is  advanced  stepwise  in  time  so  that  it  is  compatible  with 
s 

the  step  size  Ax  and  the  time  vt.  Consequently,  when  u^t  =  \x  ,  the  projec¬ 
tile  is  advanced  in  the  calculation  one  cell  dimension.  This  was  found  to 
work  satisfactorily  and  eliminated  complicated  program  logic. 

The  calculation  was  started  by  integrating  the  split  equations  in  the 
X  direction  for  a  half  time  step  using  the  initial  conditions  and  the  boundary 
conditions.  Results  from  these  calculations  were  then  used  as  initial  condi¬ 
tions  for  the  integration  for  a  half  time  step  in  the  y  direction.  The  results 
are  the  full  time  step  flow  field.  The  cycle  then  is  repeated  to  obtain  the  flow 
field  at  the  next  time  step.  The  direction  of  integration  for  the  first  and 
second  half  time  steps  may  be  interchanged  without  altering  the  validity  of 
the  procedure. 

Severe  expansion  gradients  were  encountered  along  the  axis  of  symmetry 
and  at  the  muzzle  lip  as  the  flow  turned.  In  order  to  overcome  these,  we 
found  it  necessary  not  only  to  employ  an  implicit  Godunov  method,  but  also  to 
modify  the  ADI  method  to  compute  the  flow  along  the  axis  of  symmetry.  The 
problem  along  the  axis  occurred  because  the  severe  radial  gradients  of 
independent  variables  could  not  be  accommodated  adequately  by  the  uncoupled 
X  and  y  sweeps  in  ADI.  As  a  result,  we  found  it  necessary  to  modify  the  x 
sweep  of  the  radial  velocity  (v)  equation  so  that  the  initial  radial  velocity  in 
the  cell  next  to  the  axis  was  a  linear  average  across  the  cell  rather  than  zero. 
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In  order  to  simplify  the  study,  we  neglected  the  details  of  the  flow  about  the 
projectile  nose  by  assuming  the  projectile  to  be  long  and  of  uniform  cross 
section.  Also,  the  adiabatic  gas  constant  was  assumed  to  be  constant  at  a 
value  of  y  =  1.2.  These  last  two  restrictions  can  be  removed,  but  the 
character  of  the  flow  will  not  change. 

As  mentioned  earlier,  the  flow  field  consists  of  an  outer  blast  wave,  a 

Mach  disc,  a  contact  surface,  and  jet  shocks.  It  is  difficult  to  treat  these 

phenomena  as  distinct  discontinuities,  since  the  computer  program  logic  can 

become  tedious.  Furthermore,  the  turbulent  mixings  will  smear  the  inviscid 

flow  discontinuities.  As  a  result,  the  Godunov  method  is  used  to  calculate 

12 

these  quantities  without  explicit  fitting.  Recently,  Sod  and,  previously, 

1 3 

Taylor,  et  al.  demonstrated  that  Godunov's  scheme  is  particularly  suitable 

for  computing  flows  with  embedded  compression  shocks  and  contact  surfaces. 

In  fact,  the  results  of  these  two  studies  indicate  that  the  Godunov  scheme  is 

preferable  in  most  cases  over  higher  order  finite  difference  methods  such  as 

LAX-Wendroff  and  MacCormack.  The  reason  for  this  may  be  explained  by 

14 

Moretti  in  his  recent  derivation  of  what  he  terms  the  "Lambda  scheme." 

Moretti  points  out  that  methods  which  take  into  account  the  physics  of  the  flow 
in  their  application  tend  to  perform  better  than  methods  derived  solely  on  the 
basis  of  a  mathematical  accuracy.  As  a  result,  he  suggests  that  the  Lambda 
scheme,  which  incidentally  is  very  similar  to  the  Godunov  approach  is 
superior  to  the  classical  finite  difference  methods. 

1  O 

Sod,  G.A.,  "A  Survey  of  Several  Finite  Difference  Methods  for  Systems  of 
Nonlinear  Hyperbolic  Conservation  Laws,  "  J.  Computational  Physics,  27, 
1978,  pp.  1-31. 

^^Taylor,  T.  D.,  Masson,  B.S.,  and  Ndefo,  E.,  "A  Study  of  Numerical  Methods 
for  Solving  Viscous  and  Inviscid  Flow  Problems,  "  J.  of  Comp.  Phys.,  ^(1), 
Feb.  1972. 

^"^Moretti,  G.,  "The  X-difference,  "  J.  Computers  and  Fluids,  j7{3). 

Sept.  1979,  pp.  191-206. 
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3.  COMPUTATIONAL  RESULTS 


The  mathematical  model  described  above  was  incorporated  into  a  com¬ 
puter  program  for  the  CDC  7600.  Calculations  were  performed  for  the  muzzle 
blast  from  an  M-16  rifle  for  the  purpose  of  comparison  with  the  experimental 
results  of  Schmidt,^  from  whom  the  initial  conditions  in  the  gun  were  obtained. 
Figure  4  depicts  Schmidt's  data  on  muzzle  pressure  and  flow  exit  velocity. 

For  these  data,^  the  muzzle  diameter  is  D  =  5.  56  mm  and  the  shell  velocity  is 

U  =  3100  ft/sec.  With  these  data  used  as  input,  the  inviscid  equations  were 
s 

numerically  integrated  to  obtain  the  time  history  of  the  flow.  The  time 
histories  of  the  Mach  disc  and  blast  wave  equations  are  illustrated  in  Fig.  5. 
Also  shown  are  the  data  of  Schmidt.  The  analytical  results  agree  reasonably 
well  with  the  experimental  measurements.  Note  that  the  blast  wave  motion 
appears  to  follow  a  power  law,  i.e. 

This  seems  to  support  the  strong  blast  wave  theory  which  predicts  that  the 

0  6 

spherical  blast  obeys  X  t  *  .  The  centerline  static  pressure  and  velocity 

distribution  at  t  =  230  fisec  are  given  in  Figs.  6  and  7,  respectively.  Schmidt*s 
experimental  results  also  are  shown.  A  tremendously  large  pressure  de¬ 
crease  at  the  muzzle  exit  is  discerned  in  Fig.  7.  Any  stable  numerical  scheme 
must  properly  take  into  account  this  significant  pressure  drop. 

The  pressure  profile  versus  Y/DatX/D  =  5.2  at  various  times  is 
illustrated  in  Fig.  8.  The  position  of  the  outer  blast  wave  and  inner  jet  shock 
can  be  identified  in  this  figure.  It  should  be  pointed  out  that  the  lateral  jet 
shock  location  at  t  =  352  fisec  and  233  fisec  is  very  close  (i.e.,  Y/D  ^13). 

Also,  the  numerical  results  (e.g.,  pressure)  with  a  "shock-capturing”  tech¬ 
nique  have  little  numerical  oscillations  or  wiggles  across  the  shock  and  con¬ 
tact  surface  discontinuities. 
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Fig.  6.  Velocity  Distribution  Along  the  Centerline 


PRESSURE,  psi 


Fig.  8.  Pressure  Profile  at  x/D  =  5.2 
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Figure  9  illustrates  the  velocity  vector  projections  of  the  muzzle  blast 
flow  field  at  t  =  205  fisec.  The  envelope  of  the  outer  blast  can  be  located 
easily  at  X/D  ~  31  along  the  centerline.  The  Mach  disc  and  the  associated 
triple  point  position  can  be  discerned  as  well  [e.g.,  point  ^ 

The  second  computed  results  correspond  to  a  4.2-in.  mortar.  The 
initial  conditions  inside  the  barrel  are  p  =  1300  psi,  U  =  900  ft/sec,  and 
T  =  1700°K.  The  conditions  outside  the  barrel  are  p  =  14.  7  psi  and  u=v=0. 

The  projectile  assumes  a  speed  of  900  ft/sec.  The  calculations  were  per¬ 
formed  with  and  without  a  projectile  to  determine  the  effect  on  the  muzzle 
blast.  The  pressure  and  density  contours  for  both  cases  are  shown  in 
Figs.  10  and  11.  The  results  show  that  the  existence  of  a  projectile  has  a 
significant  effect  on  the  location  of  the  blast  along  the  centerline.  However, 
its  influence  on  the  lateral  development  of  the  muzzle  blast  flow  field  is  small. 
In  Fig.  10,  the  approximate  position  of  the  projectile  can  be  seen  in  the 
density  contour. 

It  takes  approximately  12  min  on  a  CDC  7600  computer  to  integrate  the 
governing  equations  to  t  =  600  /isec  with  a  120x120  mesh  system. 
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Fig,  9.  Velocity  Vector  Projection 
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Fig.  10b.  Constant  Pressure  and  Density  Profile  for  a  4.2-in 
Mortar  with  Shell  [Pressure  (Muzzle  Blast) 

T  =  0.  203E-02] 
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Fig.  11b.  Constant  Pressure  and  Density  Profile  for  a  4.2-in 
Mortar  with  No  Shell  [Pressure  (Muzzle  Blast) 

T  =  O.I9IE-O2] 
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Fig.  11c.  Constant  Pressure  and  Density  Profile  for  a  4.2-in 
Mortar  with  No  Shell  [pxlO^  (Muzzle  Blast) 

T  =  0.  191E-02] 


4.  SUMMARY 


A  numerical  model,  which  uses  Godunov's  nonlinear  scheme,  is  formu¬ 
lated  to  solve  the  time -dependent  axisymmetric  muzzle  blast  flow  fields. 
Results  are  given  for  an  M-16  rifle  and  a  4.  2-in.  mortar.  The  calculated 
results,  including  the  effects  of  a  shell,  are  compared  with  Schmidt's  data, 
and  the  agreement  is  encouraging.  Since  the  flow  patterns  are  complex  (con¬ 
sisting  of  outer  blast  wave,  Mach  disc,  contact  surface,  jet  shock,  and  large 
expansion  fans),  our  numerical  model  with  a  simple  shock  smearing  scheme 
works  well.  This  good  agreement  is  not  considered  to  be  fortuitous.  We 
consider  that  it  results  from  using  Godunov's  scheme  which  has  properly 
taken  into  account  the  waves  propagation  and  domain  of  dependence  principle. 
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